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Abstract 

We report the results of a detailed numerical analysis of a real resonant spherical gravitational 
wave antenna operating with six resonant two-mode capacitive transducers read out by supercon- 
ducting quantum interference devices (SQUID) amplifiers. We derive a set of equations to describe 
the electro- mechanical dynamics of the detector. The model takes into account the effect of all 
the noise sources present in each transducer chain: the thermal noise associated with the me- 
chanical resonators, the thermal noise from the superconducting impedance matching transformer, 
the back-action noise and the additive current noise of the SQUID amplifier. Asymmetries in 
the detector signal-to-noise ratio and bandwidth, coming from considering the transducers not as 
point-like objects but as sensor with physically defined geometry and dimension, are also investi- 
gated. We calculate the sensitivity for an ultracryogenic, 30 ton, 2 meter in diameter, spherical 
detector with optimal and non-optimal impedance matching of the electrical read-out scheme to 
the mechanical modes. The results of the analysis is useful not only to optimize existing smaller 
mass spherical detector like MiniGrail, in Leiden, but also as a technological guideline for future 
massive detectors. Furthermore we calculate the antenna patterns when the sphere operates with 
one, three and six transducers. The sky coverage for two detectors based in The Netherlands and 
Brasil and operating in coincidence is also estimated. Finally, we describe and numerically verify 
a calibration and filtering procedure useful for diagnostic and detection purposes in analogy with 
existing resonant bar detectors. 
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I. INTRODUCTION 



Resonant bar antennas are in continuous operation at A.2K with sensitivity and band- 
width never reached before These detectors could improve their sensitivity of one 
order of magnitude in the coming years when operating at lOOmK. The first ultracryo- 
genic spherical gravitational wave detectors {jj 4] are currently completing their engineering 
phase and will soon be operational with an expected sensitivity better than 10~ 21 if z -1 / 2 
at 3kHz. Much interest is now directed towards the next generation of acoustic detec- 
tors, which will be large mass spheres [5] equipped with traditional resonant transducers 
and broadband dual resonators jf], Q]. The resonant spheres rely on available technology 



HQ, 



B S 10, 3, Q, Q, and the ground-breaking 



a rather extensive theoretical work 
experimental work performed so far on the ultracryogenic small mass sphere Minigrail 
The wide-band dual detector potentially solves the problem of relative narrow bandwidth 
of current resonant bar and spheres, but the required technology needs to be assessed in 
separate experiments [jj]. Here we numerically analyse the sensitivity and performance of an 
ultracryogenic spherical detector, equipped with capacitive, SQUID-based, resonant trans- 



ducers. A general analysis of the problem o 



given in [14J, [15( . Previous works [10 
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die read-out system of for linear detector is 
have provided a solution for the 
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read-out and inverse problem of a spherical detector, deriving equations of motion for the 
five degenerate quadrupole spheroidal modes coupled to N identical, point-like, single-mode 
resonant transducers located at arbitrary points on the sphere surface. The effects of trans- 
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ducer asymmetries on the strain sensitivity and bandwidth was studied in 
generic radial, point-like, single mode, identical resonators. In [jjj] the strain sensitivity for 
a spherical gravitational wave detector with a three-mode inductive transducer with optimal 
parameters is calculated. 

The present paper shows the results of a detailed numerical calculation of the performance 
of a spherical detector, which uses 2-mode capacitive transducers where the electrical reso- 
nant mode of a superconducting matching network can be tuned to the resonant mechanical 
modes. The signal current from the superconducting matching transformer is read-out by 
sensitive SQUID amplifiers. We chose this transduction system mainly for two reasons: first, 
most existing bar and spherical antenna use capacitive transducers, second, the technology 
involved is so far the most advanced. As recently experimentally demonstrated on the cur- 
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rently most sensitive resonant antenna AURIGA such a read-out scheme enhances both 
the sensitivity and the bandwidth of a resonant detector when working in the tuned mode. 
Two-stage SQUID amplifiers coupled to a high quality factor load can reach nowadays an 
energy resolution only an order of magnitude higher than its quantum limit when properly 
cooled down to lOOmK 20J. Two-stage SQUID amplifiers operates at 4-5K on the bar an- 



tenna AURIGA [H, and on the spherical antenna MiniGRAIL (2l|, with an unprecedently 
reached sensitivity of about 600/1. 

We use the analysis presented here to study the detector sensitivity as a function of the 
SQUIDs, the superconducting matching network and the mechanical resonators intrinsic 
parameters, and to define the optimal coupling between antenna, transducer and amplifiers. 
The simulations consider the effect of all the parameters involved in a real detector, including 
the effect of the cold damping network used in the fluxed lock-loop (FLL) of the SQUID 
amplifier to stabilize the read-out of high Q loads j^. We study also the sensitivity and 
signal bandwidth deterioration coming from the transducer being a geometrically extended 
object rather than a point-like mass. Finally we describe and numerically test a method 
to fully calibrate a spherical detector and to derived the optimal filter parameters from the 
experimental data. This method is a generalization of the one used with the resonant bar 
antenna AURIGA 



231 ] . The codes generated to perform such an analysis can be used as a 
guideline for the development of future detectors as well as a tool to evaluate the performance 
of present small mass spherical detectors. 

In section |TT] we give an overview of the equations necessary to describe the coupling of a 
gravitational wave to a sphere and of a sphere to N mechanical resonator following the for- 
malism introduced by Johnson and Merkowitz . We complete the equation of motion for 
a capacitive transducer coupled to a SQUID amplifier through a superconducting matching 
transformer. The detector strain sensitivity, noise temperature and signal bandwidth are 
derived using the generilized vector approach proposed by Stevenson 13j. This method is 
particularly powerful when transducers are not identical not only in their mechanical pa- 
rameters, but also with respect to their noise sources. It indicates a rather simple method 
to handle the correlation between transducer output channels and to form statistically in- 
dependent channels (the mode channel concept used in 16|) in absence of symmetries. In 
section II III we describe all the noise sources acting on the sphere and the transducer chain. 
The results of the numerical analysis are reported in section IIVI First, following the de- 
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scription of a parameter optimization procedure, we calculate the strain sensitivity for an 
ultracryogenic spherical detector, 2 m in diameter, made of a copper and alluminum alloy 
(CuAl6%), and equipped with one and six transducers located in the positions of the Trun- 
cated Icosahedral (TI) configuration proposed by Johnson and Merkowitz . Secondly, the 
antenna patterns are calculated for two spherical detectors with optimal sensitivity located 
in the north and south hemisphere of the earth. The sky coverage is estimated for such 
spheres operating with one, three and six resonators, both when working independently and 
in coincidence. Finally, the effect on the sensitivity is evaluated when non-identical and 
non-ideal transducers are employed. In section |V] a complete calibration procedure for a 
spherical detector is described. First we show a way to derived experimentally the equiva- 
lent temperature of the mechanical modes of a sphere equipped with six transducers. Then 
we discuss a procedure to measure the detector transfer functions and to calibrate its force 
response. Such a method is finally tested by calculating the detector mode channels response 
with simulated gravitational wave bursts. Finally, section |VT] summarizes the results. 



II. SPHERICAL GW RESONANT DETECTORS 



In this section we describe the model used to derive the sensitivity of a spherical grav- 
itational wave antenna with N transducers coupled only to radial motion. We consider a 
capacitive transducer, which can be operated in a 2-mode configuration when the electrical 
mode is tuned to the mechanical one. 

The dynamics of a sphere coupled to radial transducers and interacting with gravitational 
waves is described below using the mode channels formalism introduced by Johnson and 
Merkowitz [if]], and matrices notation suitable to derive the signal-to-noise ratio, SNR, and 
the noise temperature of a multichannel system 



A. Coupling of a bare sphere with the gravitational field 

A gravitational wave is a time- dependent deviation of the metric perturbation. In the 
coordinate frame of the wave, denoted here by primed coordinates and indices, with the 
origin at the detector center of mass and the z' axis aligned with the propagation direction 
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of the wave, the spatial metric perturbation in the TT gauge is given by 

/ h+(t) h x (t) 0\ 
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with h + and h x corresponding to the two independent wave polarizations. A wave with 
polarization h + deforms a test ring into an ellipse with axes in the x and y directions. A 
wave with polarization h x deforms the ring at a 45-degree angle to the x and y directions. 
A circularly polarized wave has h + = iL and rotates the deformation of a test ring in the 



right-handed (or left-handed) direction 241 ] . 



The dynamics of a bare resonant sphere can be described by ordinary elastic theory. A 
forces / acting on the body will cause a displacement of the sphere mass element at its 
equilibrium position. The mechanics is described by the standard equations of motion of a 
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161 ]. In particular we make use of the formalism 



Denoting by f(x, t) the total force that acts on the sphere, including the gravitational 
wave force, at the position x = {x, y, z} and time t, one finds that the equation for each 
mode amplitude is the one for a forced harmonic oscillator. After Fourier transforming, each 
mode amplitude can be written as 

a m (w) = 4p ^ ^-wa+ja^*™ x 
x J * Zm (x)f (x,w) d 3 x, 

where R is the radius of the sphere and p the density. On the right hand side of the 
equation, the first factor is an arbitrary normalization constant. The second factor describes 
the oscillating nature of the displacement where u m and $ m = 1/Q m are, respectively the 
resonance frequency and the loss angle associated with the quality factor Q m of the mth 
mode. The integral is calculated over the entire volume of the sphere where &i tm (x) are the 
time independent orthogonal elastic eigenfunctions of the sphere with I — or I — 2 . 

In general relativity only 5 quadrupolar modes of vibration (£ = 2) will strongly couple 
to the force density of a gravitational wave due to the fact that the tensor H' is traceless. 



In a perfect sphere they are all degenerate, having the same angular eigenfrequency ujq. The 
quadrupolar modes can be written in terms of the convenient set of the five real spherical 
harmonics YmiO-, d>), which are defined as follows: 
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(3) 



They are the result of a linear combination of the usual complex-valued spherical har- 
monics Y 2m . 

For a sphere of radius R the eigenfunctions can be written as: 



V m =[a{r)r + (3{r)RV]Y n , 



(4) 



The radial eigenfunctions a (r) and (3 (r) determine the motion in the radial and tangential 
directions, respectively. An explicit description of the motion in the radial and tangential 
directions is given by Ashby and Dreitlein Qj. 

In the lab-frame with origin at the center of mass of the detector and the z axis aligned 
with the local vertical, a gravitational wave produces an effective time dependent tidal force 
on each mode m of the sphere equal to the overlap integral of Eqs. (j2J). One finds 



F*(t) = ^—ph m (t)R* [cJ 2 (qR) + 3d J 2 (kR)] 
= -h m (t)m s xR, 



(5) 



where J 2 is the spherical Bessel function of order 2, the coefficients c, d specify the shape 
of the eigenfunctions and are weakly dependent on the material Poisson ratio 9j, ms is the 
physical mass of the sphere and Rx is the effective length of each mode where x depends 
on the Poisson ratio and is equal to 0.327 for the CuAl sphere considered in the following 

n 

analysis. h m are the spherical amplitudes [16J, a complete and orthogonal representation 
of the metric perturbation, q and k are respectively the longitudinal and transverse wave 
vectors as defined in The effective force on the corresponding mode of a sphere is 
therefore uniquely determined by each spherical component of the gravitational field. 
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The force acting on each spheroidal mode m in the lab frame can now be written in terms 
of the gravitational wave amplitudes 



Fm = lrn sX Rco 2 T v 




where Tv, given by 



v 



^ |(1 + cos 2 9)cos2<fi cos9sin2(f)^ 
— ~(1 + cos 2 9)sin2(f> cos9cos2(f> 



-\sin29sin$ 
^sin29cos(f> 
^-sin 2 9 



sin9cos(j) 
sin9sin<j) 




(6) 



(7) 



x 



cos2ip sin2ifj 
-sin2tp cos2%l) 



is the transformation matrix, which converts the gravitational wave amplitude in the wave 
frame into spherical amplitude in the lab frame. Here we used the y-convention of the Euler 
angles shown in figure [D and the linear combination of the spherical harmonics described in 
Eqs. ([3]). The angle if) is the first Euler angle in the rotation relating the wave frame to the 
laboratory frame and it carries information about the GW polarization. 



B. Sphere with N resonant transducers 

Resonant transducers are used on resonant detectors, either bars or spheres 
in order to improve their sensitivity and bandwidth. We consider here the same type of 
transducers as displacement sensors. They consist of a mechanical resonator with the fun- 
damental mode tuned to the quadrupolar modes of the antenna. At resonance, there is a 
transfer of momentum between the resonator and the antenna, turning small displacements 
of a large antenna into large displacements of the small resonator. 

Let us consider a set of N resonators attached to the sphere at arbitrary positions (9j, 
The values of the relative radial displacement of the sphere at the transducers location can 
be grouped together into pattern vectors for a particular mode. These column vectors, in 
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turn, may be put together to form a pattern matrix B m j defined by 26] , 

f 
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where a is the radial eigenfunction introduced in Eqs. (J4]). One gets 
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(9) 



Assuming that each resonator is designed to obey the one dimensional harmonic oscillator 
law, the coupled equations of motion for the sphere modes, written in matrix form, are 
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F s (t) 
F R (t) 



where matrices are denoted by bold fonts and capital letters and vectors by bold fonts and 
low case letters. The vector a has 5 components and the vector q has N components. They 
represent the radial displacement of the sphere and the resonator, respectively. For identical 
sphere modes and identical transducers, the factors m, k, indicating respectively mass and 
spring constants are identical and can be treated as numbers. In reality, each mode has its 
own mass, quality factor and spring constant, so they have the form of a diagonal matrix 
with components mf, kf = m^ujf(l + for the sphere, with i = 1..5 and m r - and 

kj = m T jUJj(l + i$>j(uj))) for the resonators, with j = 1..N. Here u n , with n = i,j, is the 
natural frequency and $„(co>) is the loss angle of each resonant mode. It represents the 
requency dependence of the loss of a mode. For commonly observed dissipations in metals 
27) 1. losses do not depends on frequency and $„(c^) = 1/Q n where Q n is the mode quality 
factor. In the case of viscous damping, due for example to eddy-current effects, the loss 
angle is proportional to the frequency and is given by & n (uj) = u/uj n Q n . We consider the 
first dissipation mechanism to describe the losses in the mechanical modes. F s and F R are 
the driving forces, which include the gravitational waves contribution as well as the forces 
generated by noise sources. 

Equation (fTOj) fully describes the mechanical system sphere-resonators when the pattern 
matrix B m j is known. Here we consider the special transducer configuration proposed by 
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FIG. 1: a) The truncated icosahedral (TI) arrangement for a spherical gravitational wave antenna 
with resonator locations indicated. The numbering of the resonators corresponds to the ordering 
used in the numerical calculations, b) Euler angle transforms convention. 



Johnson and Merkowitz 



lfl 



281 ] . It consists of a set of six transducers placed on the 6 



pentagonal faces of a Truncated Icosahedron (TI). The resonators are located at two polar 
angles, 6 TI = 37.3773° and 79.1876° as illustrated in figure [TJ Their azimuthal angles ipTi 
are multiples of 60°. 

Below we derive the complete equations of motion for a spherical detector equipped with 
capacitive transducer and SQUID amplifiers. In a capacitive transducer, the resonating 
mass, tuned to the spheroidal modes, modulates the charge of a parallel plate capacitor 
biased at a large constant voltage. The capacitor is formed by the resonating mass top surface 
and an electrode, assembled with a gap of the order of tens of micrometers. The input coil 
of the dc-SQUID chip is coupled to the capacitive transducer via a high-Q superconductive 
transformer, which can have, eventually, the electric resonance coupled to the mechanical 
modes in order to enhance the bandwidth. The superconducting transformer is essential 
to match the low impedance of the SQUID with the high impedance of the capacitor. In 
the model, which is schematically shown in figure [2} we include the relevant Gaussian noise 
sources of the read-out scheme. 



C. Equation of motion of a spherical detector with resonant capacitive transducers 

The motion equations of a single capacitive transducer coupled to one mechanical mode 
of the sphere can be generalized to include the complete mechanical response of a sphere 
coupled to N resonators, described by Eqs. ffTUl) . and all the equations for the electrical 
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FIG. 2: Single mode electro-mechanical model of a spherical antenna with mechanical resonator 
and capacitive transducer coupled to a SQUID through a superconducting matching transformer. 

circuit of each resonator. After Fourier transforming, Eqs. ffTUl) can be simplified denning 
the following 5 + N square matrices: 
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They are respectively the mass, elastic and force matrices of the coupled sphere. The 
force matrix A describes the mechanical coupling between the 5+N resonant modes of the 
detector. We can finally write 



[-u 2 M + K] 



a(w) 






:«)" 




= A 




q(w) 









(12) 



We assume that each of the iV transducers mounted on the sphere has the same electrical cir- 
cuit configuration described in figure [2, but not necessarily the same value of the parameters. 
We shall define then the vectors I p (k>) = (I Pi i..I Pi n), liO^) — {.U,i--U,n)^ V r (u;) = (K,i--K,jv)> 
V n (c<j) = (V^i-.V^jv), which describe, for each transducer, respectively the current in the 
superconducting matching transformer, the current in the input coil of the dc SQUID am- 
plifier, the voltage noise generated in the LC superconducting resonators and the voltage 
noise of the dc SQUID amplifiers. 

Denoting by Ei the electric field stored in the ith capacitive transducer, we define the 
electric field matrix as follows: 

E = diag{E 1 ..E N ). (13) 



10 



The electrical circuit equations can be written in matrix form as follows 



E Zn Z 12 
Z21 Z22 



I 






I 




V n (w) 



(14) 



We chose the electric matrix E to be diagonal, because we made the reasonable assump- 
tion that the electric field force acts only on the resonator which the field is applied to. The 
impedance matrix 

Zn Z 12 
Z21 Z22 

is a (N + N) x (N + N) matrix. Each of the four N xN matrices is diagonal if we consider 
as negligible the possible crosstalk between the read-out electronics of each transducer. Each 
diagonal member of Zjj is equal to 



(15) 



; 12 



-jujM i 



(16) 



Z* 21 = -juM i , Z\ 2 =juj{U s + L\ n ), 

where r % is a resistance associated with the losses in the superconducting resonator. C l is the 
total transducer capacitance resulting from the parallel between the transducer and parasitic 
capacitance. r s is a lossless resistance resulting from operating the SQUID amplifier in flux 
locked loop (FLL). There is no thermal noise contribution associated to r s because it is 
the result of a feedback mechanism. Such a resistance can be controlled by implementing a 



cold-damping system [22j. In this way the detector has a virtual low quality factor making 
the FLL electrically stable. We remark that introducing such a damping scheme in our 
calculations brought benefit to the numerical analysis, by eliminating the computational 
problem of sharp resonances. 

To fully describe the mechanical and electrical dynamics of the detector we have to 
introduce the back-action of the electrical read-out circuit on the mechanical system. The 
current flowing in the LC loop of the circuit in figure [2] generates a force on the mechanical 
resonator proportional to the current itself and the applied electric field. The back-action 
force vector Fjj A (u;) adds to the Langevin force generators Fj$(o;) introduced in Eqs. ( |T2l) 
and is equal to 
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The complete set of coupled equations of motion becomes finally: 



M Z BA 
E 

Z 





where M. = — cj 2 M + K, Zba is the (5 + N) x N back-action matrix given by 
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J BA 
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and 



A' 



A 
I 



(20) 



The (5 + 3iV) square matrix on the left side of Eqs. (|18[) can be seen as the impedance 
matrix Z of the electro-mechanical system. Defining G = Z~ l A', the SQUID input current 
for each transducer is given by 



Ii = G, 



f£M 



GiF, 



(21) 



V r (w) 

where Gi is a submatrix of the admittance matrix G with components G Sjr , where s = 
5 + 2N..5 + 3N and r = 1..5 + 3iV. The vectors F s , F R ,V r and V n are the forces generated 
by each noise source. 

The noise of the detector, referred to the SQUID amplifier input, in absence of signal, is 



described by the spectral density matrix Si 



29j . Each component of the matrix, 



(22) 



is the Fourier transform of the correlation function for the mth and nth outputs defined as 

/"OO 

Rnm(r) = (IT(t)I?(t ~r))= If{W(t - r)dt. (23) 



From Eqs. (123]) . each component of the spectral density matrix becomes the product of the 
Fourier transforms of the transducer outputs, as shown in the second equality of Eqs. (T221 . 
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The white current noise 1$ of the SQUID amplifier, which will be better defined below, 
needs to be added to the spectral matrix Si. In matrix notation we can easily write the 
total SQUID current spectral density matrix as 

Si = GiFF*Gj + S Ij0 , (24) 

where the N square matrix Si.o has components S™^ 1 = I™Iq*. The diagonal elements of 
SV,o are equal to the current spectral density given in Eqs. ( 1361) . The correlation between 
the iV SQUIDs additive current noise is expected to be negligible so, in the following, 
we will consider diagonal the matrix Sio- The spectral density defined in Eqs. (124"]) can 
be numerically calculated and experimentally measured by means of techniques where the 
phase information is preserved. 

The optimal signal to noise ratio po , f° r a gravitational wave signal of amplitude h(ui), 
is given by 

J S hh {uJ) Ztt 

where 

S hh {u) = h» [(F^)*G: igiI S I - 1 (^)G sig , I F^]' 1 h(^) 

= {m s xR^f [ T vG: ig ,iS I -(-)G sig , I T v ]- 1 (26) 

is the one-sided total strain noise power spectrum. 

In the equation above, we called h the vector (h + h x ) and the N x 5 matrix G s i g> i is a sub- 
matrix of the admittance matrix G(u;) with components G s ^ r , where s = 5 + 2N + 1..5 + 3N 
and r = 1..5. 

Each transducer line should be considered as a linear system with the (5 + 2A^ + 2) uncor- 
rected noise sources described above, if we consider as negligible the correlation between the 
voltage and current noise in the SQUID amplifier. However, the outputs of the N transducers 
do have correlated noise and the off-diagonal components of the spectral density matrix are 
non zero. One can always find linear combinations of transducers outputs, which produce 



N uncorrelated signals [13|. Since Si is an N x N Hermitian matrix it can be diagonalized 
by an unitary matrix U(oS). A new output channel vector is then obtained and is related to 
the original vector Ii by 

II 1 = U*!;. (27) 
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The channels I" are statistically independent and the spectral density matrix 

= ITSiU = diag{il) X {w), ^H) (28) 

is diagonal, the eigenvalues ^ of S" being the noise spectral density of each independent 
channel. 

We notice that, after performing the diagonalization in Eqs. (1281) . the total optimal SNR 
can be written as the sum of the SNR in each statistically independent output channel. 
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From Eqs. (|25|) and Eqs. (|26|) the SNR becomes 
'm s xRool\ 2 |M^)| 2 
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(30) 



The performance of resonant detectors is often characterized by their sensitivity to im- 
pulsive burst signals which vary little over the detection bandwidth. For impulsive signals, 
;he SNR is proportional to the deposited energy E in the antenna initially at rest given by, 



E 



7T PsVg 



U\h(u) 



(31) 



G167T Ul ^ u/l 4 f 
where S = ^^fH is the integrated cross section of a spherical detector [9j and II is the 
reduced energy cross section equal to 0.215 for a CuAl sphere [jjj]. G and c are Newton's 
gravitational constant and the speed of light respectively, while ps and vs are the antenna 
material density and sound velocity. 

The pulse detection noise temperature T N is then defined as 

T N = ±— 2 . (32) 

This is a convenient quantity to compare spherical detectors with bar detectors: while 
E for a bar-antenna depends on source direction and polarization, T/v does not. Using 
Eqs. ( !30l and Eqs. ( !3T|) we can write T/v as follows 



T, 



4ir 
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4 / TVG: igiI (^)S I - 1 (^)G sig> i(^)T v ^ 
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■ (33) 
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We remark here that, due to the dependency of the matrix Ty from the wave direction 
and polarization, the sensitivity of a spherical detector will be isotropic over the sky only if 
a sufficient number of transducer is used (N > 5). For N < 5, one can define the detector 



sensitivity by averaging over the direction and polarization as described in 



30|. 



III. NOISE CONTRIBUTIONS 
A. Mechanical resonators 

Thermal noise is the main contribution of the mechanical resonators to the total detector 
noise. The spectral density of the thermally activated forces acting on the mechanical modes 



can be estimated from the fluctuation-dissipation theorem [3l| , and described as follows 



S F , th (u) = 4k B TRe(j^,Y (34) 

where 4£ is the system mechanical impedance matrix derived in Eqs. ffl8|) . kg is the Boltz- 
mann constant and T is the thermodynamic temperature. With this formalism we take into 
account the stochastic crosscoupling among the various degrees of freedom of a macroscopic 



mechanical body 



32 ] . 



The thermal noise contribution of other modes of the spherical detector has never been 
considered as a possible source of noise because they lay generally far a way from the detec- 
tor bandwidth. In a real detector, the first toroidal modes are only a few tens of Hz lower 
than the spheroidal modes due to the spread of the resonances caused by the detector asym- 



metry 33]. Moreover, in a spherical detector with large bandwidth, the toroidal modes will 
unrecoverably fall into the sensitive bandwidth. However they are not sensitive to GW and 
they couple very weakly to a resonator with radial sensitivity. The toroidal modes usually 
have large quality factors than the spheroidal modes and their thermal noise contribution is 
generally negligible. 

Further-more, higher frequency modes could also contribute to the total noise in the de- 
tector bandwidth due to down-conversions phenomena related to the physical geometry and 



dimension of the read-out transducer 



34| . To evaluate the contribution of higher frequency 



mode one can proceed as in 35j , and calculate the total mechanical impedance using Finite 



Element or numerical techniques. Here we consider only the thermal noise of the spheroidal 
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modes and the main radial resonance of the transducers, the latter being usually the dom- 
inant source of noise. The contribution from other modes will be studied in a following 
paper. 



B. Electrical resonators 

The LC resonator which derives from transducer capacitance and the primary coil of the 
matching transformer contributes to the total noise with a thermal voltage noise source asso- 
ciated with the resonator losses, with single-sided spectral density Sy = Ak,BTRe{Z LC {uj))- 
The dissipating term r = Re{Z^c{^)) is linked to the intrinsic electrical quality factor of the 
LC resonator by the well-known relation Q e i = u res L p / 'Re(Z lc) • The dissipation resistance 



r includes the contributions from dielectric losses in the transducer and decoup 
in the coil parasitic capacitance and in the coil insulating layers and holder 



ing capacitor, 
magnetic 



losses due to flux motion in the superconducting s 



lields, "magneto- resistive" losses due to 



dissipative components in the SQUID amplifier 37]. 

The thermal noise contribution from the LC resonator adds to the SQUID back-action 
noise and may become significant when SQUID amplifiers with e < 200ft are employed. 

C. SQUID amplifier noise theory 

A coupled dc SQUID amplifier can be modelled as an ideal current amplifier with a 
current noise source I n in parallel and a voltage noise source V n in series with the input coil. 
The two of them are responsible respectively for additive and back action noise. To fully 
characterise the SQUID it is necessary to estimate both the noise contributions. An useful 
parameter to characterize an amplifier is its noise temperature T n , defined by 

rp _ \/S vv Sii - Im(S iv ) 2 _ u I ~ 

in the classical limit when /c_bT/v >> hoi. T n is the temperature at which the optimal input 
impedance gives the thermal noise power equal to the amplifier noise. Its minimum value 
'or a linear amplifier is imposed by the uncertainty principle and is given by T^simhw jks 
. S vv and Sa are the spectral densities, referred to the SQUID input coil, of the two 



noise generators, while S V i is the cross correlation between the two. According to Clarke- 



1(3 



Tesche-Giffard (CTG) theory the single-sided spectra are equal to 

11 R s \ M ilS Q J 

S vv ~ IiML^sq? (36) 
S jv ~ U^juLsq, 

where i? s , Lsq, M it sQ are respectively shunt resistors, self inductance and input coil mutual 
inductance of the SQUID. The cross correlation power spectrum is always purely imaginary 
due to the time-reversal simmetry of the SQUID motion equation. Its absolute value is 
usually small and we will neglect it in the following calculations. In Eqs. fl35l) the noise 
temperature is also defined in terms the energy resolutions en, e vv and e„ expressed in units 
of h. In the following sections we will not make distinctions between the voltage and current 
energy resolution, e„ and e vv , and we will generically talk about SQUID energy resolution 
assuming that they both have the same value. The cross correlation energy resolution e iv 
will be neglected. 



IV. NUMERICAL ANALYSIS 



Here we numerically calculate the noise temperature, the bandwidth and the strain sen- 
sitivity of an ultracryogenic, large 2m in diameter, spherical detector as a function of the 
transduction chain parameters. We consider a detector operating according to the state-of- 
art resonant antenna technology and when all the parameters are improved to operate the 
detector nearly at the quantum limit. Further we calculate the antenna patterns and the 
sky coverage of two identical detectors operated in coincidence and located respectively in 
Leiden (The Netherlands) and Sao Paolo (Brasil), the location of the two small spherical 
antenna Minigrail and Mario Schenberg 4], currently under development. Finally the 
anisotropy in the sensitivity and bandwidth is studied for a not ideal resonator and a not 
optimally tuned 2-mode capacitive transducer. 
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A. Parameters optimization 



To optimize a resonant detector one must know the voltage and current noise of the 
available SQUID amplifier. The other parameters can then be adjusted to achieve the 
best SNR. The optimal impedance matching between the mechanical resonators and the 
SQUID amplifier is achieved when the transducer electrical mode is tuned to the mechanical 
resonances. In the scheme described in figure El the electrical mode contributes to the 
transducer chain with an equivalent mass m eq = CtE 2 /uj 2 ^ In a multimode detector the 
energy transfer between each resonator is optimized when /i = mji/m e ff = m e i/m^ where 
fi is defined as the mass ratio and m e ff is the effective mass ratio of the five spheroidal 
modes, which is a fraction of the sphere total mass m s and is equal to m e ff = 5/6xm s 
For optimal mass-ratio one gets a shorter energy transfer time between the modes 
and a wider system bandwidth. To fully describe the sensitivity of a resonant detector we 
consider the signal-to- noise ratio SNR, the pulse detection noise temperature Tjy and the 



signal bandwidth, here described by ; 



39) 



Sf _ i (/o°°l^V)M 2 m 



In the analysis below we imposed the following conditions: 

1/2 



E = IMJr 



C T 



(38) 

2 



U el — C T L r ~ U tn 

where L r = L p (l — k 2 L s /(L s + Li)) is the reduced inductance of the primary coil of the 
matching transformer with coupling k 2 = M 2 /(L s L p ). The first condition is obtained by 
assuming a constant mass ratio fj, = m e i/rriR. With the second condition we impose that 
the electrical resonances are tuned to the transducers mechanical resonances. This condi- 
tion is fundamental to get optimal matching in a capacitive transducer. The parameters of 
the electrical matching network depend on the value of the SQUID input inductance and 
noise, and on the detector mass ratio. It can be shown that for large mass ratio, a cou- 
pling k ~ 1, and optimal matching, the inductance L s of the superconducting transformer 
secondary coil should be as large as the SQUID input inductance Lj. However in practice 
large L s will guarantee large coupling and finally a better sensitivity. For low mass ratio, 
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the impedance seen by the transformer primary coil is lower than for high mass ratio when 
the same capacitance is considered. One can obtain a better matching by decreasing both 
the inductances L p and L s of the primary and secondary coil while maintaining valid the 
conditions Eqs. (1351) . The improvement is however rather limited giving for example about 
20% better noise temperature using an inductance L s ~ 0.3Lj for a mass ratio fi = 0.001 
and a k ~ 1. 

We anticipate that for a detector with significant thermal noise (high T/Q ratio) oper- 
ating with a low noise SQUID amplifier, the optimal sensitivity is obtained for large mass 
ratio fi and, consequently, high transducer capacitance (see Eqs. ( |38|) ). It is however physi- 
cally impossible to reach arbitrary resonator high mass and capacitance using the detector 
configuration discussed here. The maximum transducer capacitance Ct and mass m# one 
can reasonably obtain, which however have never been developed so far, are estimated to 
be CT,max — 30mF and m^ max ~ 90Kg. They correspond to a transducer with an area 
A ~ 0.1m 2 , a gap d ~ 20/im and a mass ratio fi ~ 0.01[46(. We notice that strong physical 
constrains apply in general also to the mechanical resonant frequency due to the limits on 



the resonator membrane thickness 



40 



411 ] . However, a resonator "rosette" design as devel- 



oped for the NAUTILUS and EXPLORER detectors [4l|], allows a rather wide freedom in 
the choice of the resonator mass and sensitive area for a give resonant frequency. 

We report below the calculated detector sensitivity for a 2m large in diameter, 30 ton, 
CuAl sphere. The main parameters of the detector are summarized in table [J Such an 
antenna has a cross section X = 9.76 • \0~ M m 2 Hz. The total energy deposited by a GW 
of amplitude h(u) is E = 1.0 • I0 35 u) 2 \h(uj)\ 2 K. A GW burst signal, lasting for a time 
tq lms shorter than the detector integration time and rising quickly to an amplitude ho = 
h(u)o)/TG ~ 2 • 10~ 21 , deposits an energy Tqb = lfiK. 

We consider a sphere equipped with six radial capacitive resonators in the TI config- 
uration. The electrical mode of the superconducting matching network is tuned to the 
mechanical ones. The signal of each transducer is amplified by a two-stage SQUID ampli- 
fier. Two situations have been considered. In the first one the detector parameters have 
values according to the current available technology. In the second we estimate the ulti- 
mate sensitivity of the detector obtainable with quantum limited amplifiers and ultra high 
Q mechanical and electrical resonators. For a given T/Q ratio and SQUID amplifier energy 
resolution, the detector sensitivity depends on the resonators mass ratio, the electric field 
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current 


QUCLTltUTTl 


parameters 


J. U 1 

technology 


limit 


sphere mass , M s [tons] 


30 


30 


sphere diameter [m] 


2 


2 


resonator mass, M r [Kg] 


10 


90 


spheroidal modes, fi [Hz] 


987, 1001, 1008, 


987, 1001, 1008 




1012,1017 


1012, 1017 


QsiQriQel 


2 • 10 6 


5 • 10 7 


T [mK] 


50 


20 


C t [nF] 


10 


30 


L p [H] 


3.6 


1.2 


bias field, E [Volt/m] 


5 • 10 6 


4- 10 7 


SQUID sensitivity, E res [h] 


50 


1 



TABLE I: Main parameters used in the numerical analysis of the sensitivity of a spherical detector 
equipped with six capacitive transducers. The first column shows the parameters of a 2 m sphere 
in CuAl already achieved in separate experiments. In the second column we give the parameters 
necessary to operate the detector at the quantum limit. The CuAl alloys have a sound velocity 
of 4700m/s and a Poisson ration of about 0.3. The material dependent factor x introduced in 
Eqs. (0)is x = 0.327 for CuAl. The effective mass of the spheroidal modes is then m e ff = 
5/6x * m s ~ 8 tons. 



bias E, the transducers capacitance and the parameters of the matching transformers. 

The detector effective noise temperature, the bandwidth and the SNR for a GW burst of 
Tgb = 10fiK are calculated as a function of the mass ratio for different values of the T/Q 
ratio. T/Q refers both to the mechanical and electrical resonances. The equivalent mass 
of the electrical resonator is adjusted according to the mass ratio by changing the electric 
field bias E. The optimal electrical field is calculated as a function of the mass ratio and the 
transducer capacitance. The results are shown in figure [31 We considered mechanical and 
electrical T/Q of 1 • 10" 7 , 2.5 • 10" 8 , 4 • 10" 10 . Mechanical T/Q of the order of 10" 8 have been 
reached in a CuAl sphere cooled at 50mK 42|, and in A15056 bars cooled at 100 mK. The 
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FIG. 3: Optimal bias electric field for different value of the transducer capacitance as a function 
of the resonators mass ratio \i. 

lowest T/Q = 3- 1CT 8 for electrical resonators has been achieved by the AURIGA group with 



a large Nb coil resonator cooled down to 50mK 20| . The fabrication of electrical resonators 
with T/Q = 1 ■ 10~ 7 at acoustic frequencies is nowadays a well established technology. A 
T/Q < 10 -8 has not yet been achieved experimentally either in mechanical transducers or in 
electrical resonators. The results of the calculations are shown in plots figure Hb.) , where the 
noise temperature, the SNR and the bandwidth are given as functions of the mass and the 
T/Q ratio. The SQUID amplifier coupled energy resolution was chosen to be e coup i e d = 50/1. 
A coupled energy resolution as low as 27 h has been recently obtained with a two-stage 



SQUID amplifier coupled with an electrical resonator and cooled down to T = 50mK [20 ] . 
The spherical GW antenna Minigrail and the AURIGA detector are currently operating at 
5 K with SQUID amplifier energy resolution of the order of 600/1 H]- More than a factor 
of 10 improvement is expected when this detector will operate at T < lOOmK. 

From the plots in figure H] we see that, in terms of noise temperature and bandwidth, 
a large mass spherical detector developed using the available technology will not perform 
better than an ultra-cryogenic bar detector operating at the same frequency. As a matter 
of fact, the noise temperature T/v and the bandwidths depend only on the frequency, the 
electro-mechanical impedance matching and the SQUID amplifier noise. However, due to 
the larger cross section of a spherical detector, with respect to a bar detector at the same 



frequency, the spherical detector improves the SNR by a factor of about 40 12] . Moreover, 
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FIG. 4: Noise temperature, SNR and bandwidth of the spherical detector as a function of the 
resonator mass ratio \i. The detector is read-out by a 50h energy resolution (a) and quantum 
limited (6) SQUID amplifier. Three different configurations with T/Q = 1 ■ 10~ 7 K, 2.5 • 10~ 8 K 
and 4 • 10~ 10 K are studied. The SNR is calculated for a Tqb = lOfiK gravitational wave burst. 

omnidirectionality will of course still be a unique feature of a spherical detector when at least 
six resonators are used. With T/Q = 2.5- 10~ 8 K and mass ratio \x = 0.001, corresponding to 
a transducer resonating mass of about 9 Kg, a capacitance of 10nF and an optimal electrical 
field E = 5- 10 6 V/m, the detector has a noise temperature T N of about IO/iK, about a factor 
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of 30 better than the present most sensitive resonant bar detector 

The SQUID amplifier is considered strongly coupled to the electrical matching circuit 
assuming k = 0.8, L p ~ 1H and L s ~ Lj = \.l\iH. As in bar detectors, the sensitivity will 
be mainly limited by thermal noise of the electrical resonator, and, outside the resonances, 
by the SQUID amplifier additive noise. The curves in figure HJa clearly show that large 
mass ratio brings benefits only to the bandwidth and not to the detector noise temperature. 
This is due to the fact that a large mass resonator will not help to decrease the thermal 
noise contribution. As a matter of fact, the contribution of the mass at the denominator of 
the Langevian forces derived in Eqs. ( 1341) is almost cancelled by the transfer function when 
transforming the forces into displacement. On the other end, the bandwidth will increase for 
large mass ratio due to a better impedance matching. The latter will unavoidable enlarge 
the contribution of the back action noise coming from the SQUID amplifier producing the 
observed saturation in the detector sensitivity at large mass ratio. 

To fully exploit the potentiality of a massive spherical detector one needs to develop ca- 
pacitive transducers with high sensitive area, massive mechanical resonators with Q-factors 
larger than 10 , high Q electrical transformers and large bias electric fields. 

In figure Hlb the detector sensitivity is given for the antenna and transducer chain 
parameters obtainable by pushing the current technology to its limits. We consider a 
quantum limited SQUID amplifier operating at temperature T < lOOmK. Again we 
estimate the noise temperature T/v, the bandwidth and the SNR for a lOfiK burst for 
T/Q = 1 • 10~ 7 K, 2.5- 10 _8 i^, 4- lO -10 ^. By comparing those plots with the ones previously 
discussed, it becomes clear that without an improvement of the mechanical and electrical 
resonators with respect to available technology, the use of a quantum limited SQUID ampli- 
fier will not benefit the detector sensitivity. By looking at the dot-dashed curves of figure HI 
for examples, corresponding to a T/Q = 1 • 10 _7 i^, one shall expect no improvement in 
the strain sensitivity, but only a small increase of the bandwidth due to the lower addi- 
tive noise of the SQUID. A spherical detector operating at T = 20 mK with mechanical 
and electrical quality factor Q ~ 5 • 10 7 , a quantum limited SQUID amplifier and large 
mass mechanical transducers, \i = 0.01 and ttir = 90 Kg, can have a noise temperature 
of 1.3 • lO -7 ^, corresponding to a peak strain sensitivity of lO^^Hz" 1 ^ 2 at IkHz, and a 
bandwidth of about 200Hz. This would improve of a factor 50 the sensitivity of an existing 
bar or small sphere antenna working at the quantum limit. The bandwidth will have only 
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a moderate increase with respect to present resonant bar antennae [lj], merely due to the 
spreading of the spheroidal modes. The minimum achievable antenna noise temperature is 



given by the quantum- mechanical limit of a linear motion detection derived by Giffard [43] 
and equal to TJv TO m = 2Huj j kB[exp(Huj j ksT^^) — 1]> where T^ >a is the noise temperature of 
the linear amplifier. When a quantum limited amplifier is used, T^^ min = 9.5 • 10~ 8 K for a 
kHz resonant sphere considered here. 

B. Strain sensitivity 

We show here the strain sensitivity calculated using Eqs. fl26|) . The strain curves are de- 
rived for each noise contribution described in section UTTl The read-out circuit was optimized 
as discussed above. 

In figure [5] the strain sensitivity is calculated for a spherical detector equipped with a single 
transducer placed in the position 1. It was calculated for an optimally oriented source, as 
will be described in detail below. The detector parameters used for this simulation have 
already been achieved in separate experiments ma^, . The sensitivity curves are obtained 
considering the detector operating at T = 50mA' with a 50ft SQUID amplifier and with 
the electrical mode coupled to the mechanical ones. The quality factors are chosen to be 
Q = 3 x 10 6 and Q = 2 x 10 6 for the mechanical and electrical modes respectively. The mass 
ratio is /i = 0.001, corresponding to a mechanical resonator mass tur ~ QKg. The optimal 
electrical bias field is about E = 5 • 10 6 V/m and the transducer capacitance is lOnF. 

The strain sensitivity shown in figure [6] was calculated for a spherical detector with 
six transducers in the TI arrangement. The detector and transducers parameter are the 
same used in the single transducer configuration described above. Differently from the 
single transducer configuration, the sensitivity is independent of the wave direction and 
polarization. In such a configuration the detector sensitivity will be limited by the thermal 
noise contribution from the mechanical and electrical resonators and by the back action 
noise of the SQUID amplifier due to the the optimal matching coming from the tuning of 
the electrical resonances with the mechanical ones. Outside the resonances, the sensitivity 
is limited by the additive current noise of the SQUID amplifier. 

In figured the strain sensitivity is shown for an optimized detector equipped with SQUID 
amplifiers operating at the quantum limit. The simulation considers a detector operating 
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frequency [Hz] 

FIG. 5: Strain sensitivity for a spherical detector with one single transducer, T/Q ~ 2.5 x 10~ 8 , 
resonators mass ratio fi = 0.001 and 50h energy resolution SQUID amplifiers. The transducer 
electrical mode is tuned to the mechanical ones. The parameters used for this simulation have 
already been achieved in separate experiments with bar detectors or a lower mass spherical antenna. 
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FIG. 6: Strain sensitivity for a spherical detector with six transducers, T/Q ~ 2.5 x 10~ 8 , resonators 
mass ratio \i = 0.001 and 50h energy resolution SQUID amplifiers. The electrical modes are tuned 
to the mechanical ones. The parameters used for this simulation have already been achieved in 
separate experiments with bar detectors or lower mass spherical antenna. 

at T = 20m.fr with all resonators having quality factors Q = 5 x 10 7 . We consider six 
transducers in the TI arrangement. The mass ratio is fi = 0.01 and the optimal electrical 
bias field is E ~ 3 ■ 10 7 V/m. We consider a transducer with a large capacitance of about 
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frequency [Hz] 

FIG. 7: Strain sensitivity for a spherical detector with T/Q = 4 x 10 -10 , resonators mass ratio 
fj, = 0.01 and lh energy resolution SQUID amplifiers. The electrical modes are tuned to the 
mechanical ones. 

C = 30nF. The sensitivity is mainly limited by additive and back action noise of the SQUID 
linear amplifiers, whose minimum energy resolution is imposed by quantum mechanics. De- 
spite the low working temperature and the high Q considered, a small contribution from the 
thermal noise is still present. 



C. Antenna pattern 

We calculate here the SNR of a spherical detector equipped with one, three and six 
resonators as a function of the direction and polarization of the incident wave. We consider 
a detector operating at nearly the quantum limit. We numerically calculate the SNR as 
a function of the three Euler angles (ip, 6, <fr) of Eqs. ([7]). All the three angles are necessary 
to completely define the GW source. The first Euler angle ip carries information about the 
wave polarization 4^]. The other two, (6*,0), give the source direction. In the sky maps 
presented in this section, we call those angles respectively declination 5 = n/2 — 9 and right 
ascension a = to match the astronomical notation. The antenna patterns depend on the 
longitude and the latitude of the detector location as well as on the Universal Time (UT), 
r, due to the earth proper rotational motion. To simplify the discussion we consider here 
detectors at r = 0. We calculate the antenna patterns for two spherical detectors, 2m in 
diameter, located at Leiden (The Netherlands), Ion = 4° 30", lat = 52° 7", and at Sao Paulo 
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a) -90° b) -90° 

FIG. 8: Antenna pattern for spherical detectors respectively located in Leiden a) and Sao Paulo 
b) with a single transducer on position 1 of the six positions of the TI arrangement. We consider 
a linearly polarized (top) and a one-cycle, circularly polarized (bottom) signal as explained in the 
text, depositing an energy Tqb = 3/iK in the detector. 




FIG. 9: Antenna pattern of spherical detectors respectively located in Leiden a) and Sao Paulo 
b) with a three transducers set in the TI configuration at Oti = 37.3773° (sky maps above) and 
9ti = 78.1876° (sky maps below). We consider a one-cycle circularly polarized sinusoidal signal 
depositing an energy Tqb = 3^ K in the detector 

(Brasil), Ion = -46° 38", lot = 23° 34", where respectively the MiniGRAIL and the Mario 
Schemberg, 65 cm large in diameter, spherical detectors are being developed. We chose 
those locations because the general discussion about direction sensitivity is independent on 
the detector size and mass a part from the absolute value of SNR and bandwidth. The 
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analysis presented here is then useful for the existing small spherical detectors as well. 

In figure [8] the sky maps are shown for a spherical detector with a single transducer on 
position 1 of the TI arrangement. We consider a linearly polarized wave with h + = h and 
h x = and a random polarization angle ip (skymaps at top). The amplitude ho is related 
to the deposited energy E = k^TcB according to Eqs. (13"Tj) . For each simulation we indicate 
the GW energy in the figure captions. The SNR for each wave direction depends on the 
polarization angle. The choice of a random polarization for each direction trial produces 
the scattered pattern. When a circularly polarized wave is chosen (skymaps at bottom) the 
maps become smoother and the SNR is only dependent on the wave direction. Here and 
in the following we consider a one-cycle, circularly polarized sinusoid: h + = h,Q\/2cos(u)t) 
and h x = h(,y/2sin(ujt) , < tut < 2ti, with the frequency uj laying within the detector 
bandwidth. Transducers on other positions will show a similar pattern rotated of proper 
angles accordingly to the positions. The sky maps on the left (figure [HJ a) refers to a detector 
in Leiden with the lab-frame oriented so that the z-axis is pointing to the local vertical and 
the x axis to the local south. The sky maps on the right (figure Elb) refers to the a detector 
in Sao Paulo. As expected, in a spherical detector operating with only one transducer the 
sensitivity is direction and polarization dependent and changes according to the transducer 
location. A source emitting a linearly or circularly polarized wave is more likely to be 
detected when laying on the plane perpendicular to the transducer axis and passing through 
the sphere center. 

In figure [91 the sky maps are shown for a spherical detector with three transducers placed at 
the three positions of the TI arrangement, respectively at 6ti = 37.3773° and 9ti = 78.1876°. 
The signal used for the simulation is a one-cycle circularly polarized sinusoid with energy 
Tqb — 3 fiK. The higher maximum SNR with respect to the single transducer configuration 
is mainly due to the large detector bandwidth. When six identical resonators are used, the 
detector becomes independent form the wave incoming direction and polarization 44j . This 
result is also shown in figure fTUl For a real detector with identical transducers, not optimally 
tuned to the 5 spheroidal modes, a small residual anisotropy in the SNR and bandwidth 
is still present. This is however less than 10% over the whole sky. The signal bandwidth 
anisotropy is shown in figure [11] Both the figures will be further discussed in the following 
section. 

Figure fT2l presents the detection efficiency of a spherical detector equipped with only one 
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FIG. 10: SNR of a quantum limited spherical detector with six transducers in the TI configuration. 
In a), identical transducers are considered. In b), we arbitrary modified the parameters of each 
transducer of a maximum of 15% from their optimal value. We consider a one-cycle circularly 
polarized wave with energy Tqb = 3fiK. 
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FIG. 11: Signal bandwidth of a quantum limited spherical detector with six transducers in the 
TI configuration. In a), identical transducers are considered. In b), we arbitrary modified the 
parameters of each transducer of a maximum of 15% from their optimal value. We consider a 
one-cycle circularly polarized sinusoidal signal with energy Tqb = 3[iK. 

transducer as a function of the GW incoming direction for different GW energies, Tqb = 
"., 1.5, 3 \iK equal to a maximum SNR respectively of p max ~ 10, 15, 30. In analogy with 
451 ] . the detection efficiency is defined as ^erfc[(r) — po)/v / 2], where i] is a threshold chosen 
equal to 5 and erfc is the complementary error function. The detection of a one-cycle 
circularly polarized wave of energy Tqb = 3 fiK is almost always succesfull independently 
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FIG. 12: Comparison of detection efficiency sky maps for a spherical detector with one transducer 
and for different deposited energy Tqb = 1, 1-5, 3 fiK (from top to bottom). We consider a 
circularly polarized wave as described in the text. 




FIG. 13: Comparison of detection efficiency sky maps for a spherical detector with three transducers 
in the TI positions respectively at Oti = 37.3773° and 6ti = 73.1876°, and for a complete detector 
with 6 transducers (from top to bottom). We consider a GW burst deposited energy of 0.6 [iK 
and a circularly polarized wave. Note the difference in the color code on the graphs and that the 
signal amplitude is smaller than the ones in figure [T2l 

of the incoming direction. In figure [TBI similar skymaps are shown for a Tqb — 0.6 [iK when 
a spherical detector with 3 and 6 transducers is considered. The results of the skymaps in 
figure [12] are summarized in figure HH where the fraction of sky is plotted as a function of 
the detection efficiency for a single transducer configuration and different signal amplitudes. 
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FIG. 14: Fraction of sky as a function of detection efficiency for a detector with a single 
transducer. We consider a circularly polarized wave with Tqb = 0.6, 0.75, 1, 1.5, 3pK, and 
Pmax ~ 6, 7.5, 10, 15, 30, respectively denoted by a, b, c, d and e. 

One notes that the fraction of sky covered decreases when the detection probability level 
increases: the curves evolution for each SNR can be understood from the patterns of the 
skymaps in figure [12j The detection probability is higher than 30% in 40% of the sky for 
a Pmax = 10, corresponding to a 1ms GW burst of amplitude h = 2 ■ 10~ 21 ; it is 50% for 
almost a 70% fraction of the sky when p max = 15. For p max ~ 30 the detection is likely 
almost in any direction. 

Figure [15] presents the fraction of sky as a function of the detection efficiency of a GW 
burst of deposited energy Tqb — 0.6piK for different transducer configurations. Plot i) 
refeers to a single spherical antenna, while plot ii) shows the detection efficiency of two 
optimally oriented spherical detectors located at Leiden and Sao Paulo. This was obtained 
by maximizing the portion of sky simultaneously seen by both detectors as a function of 
the orientation of the single detector reference system with respect to the local south. This 
corresponds to a rotation of the brasilian detector lab frame (x, y, z) (see figure [I]) around 
the local z-axis of about 0o = —135°. 

A detector with a single trasducer is unable to detect any signal from the sky when 
the SNR is as low as p max = 6. For the same SNR, the sky coverage of a single detector 
with three and six transducer is respectively more than 40% and 100% with 50% detection 
efficiency. This is mainly due to the increase of detector bandwidth when more than one 
transducers are used. Figure [TSjii shows that for detector with one or three transducers, 
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a) single detector 



b) two detectors 




20 40 60 80 100 20 40 60 80 100 

detection efficiency [%] detection efficiency [%] 



FIG. 15: Fraction of the sky as a function of the detection efficiency for a GW burst of deposited 
energy Tqb = and p max = 6 for a single spherical antenna i) and two optimally oriented 

spherical detectors ii) located at Leiden and Sao Paulo. Curve a corresponds to a read-out config- 
uration with a single transducer located in position 1 of each detector, curves b and c to a three 
transducers configuration with Oti = 37.3773° and 6ti = 73.1876° and curve d to a complete 
spherical detector with 6 transducers. 

the twofold coincidence probability remains lower than the detection efficiency of a single 
sphere. Two spherical antenna with six transducers can equally detect any sources in the 
sky even when the SNR is as low as p max = 6. 



D. Sensitivity of a sphere with not ideal transducers 

It has been shown that for a perfect sphere the resonators mistuning and misplacing has 
little effect on the isotropy when the deviation from the ideal TI configuration is less than 
1% 13, . We study here two other possible degradation effects which may arise in a real 
spherical detector. The first is related to the broken spheroidal mode degeneracy and the 
second to the the fact that a real transducer is not a pointlike mass as generally considered. 
In a real detector the spheroidal mode degeneracy is broken due to the suspension and the 
holes made on its surface to house the transducers. The modes spreading can be as large 
as 5% of the main resonance. It becomes natural to ask to which modes each transducer 
should be tuned and how good the tuning should be to avoid a sensitivity degradation. We 
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consider here a numerical analysis of the transducer mechanical and electrical mistuning for 
a sphere with six resonators with a mass ratio /j, = 0.01 and operating at nearly the quantum 
limit. 

In figure [TBa) . we see the effect on the SNR for a Tqb — 3fiK circularly polarized GW 
burst when the natural resonances of each resonator are modified from the initial values 
ujq = 27r[987, 1001, 1001, 1008, 1012, 1017], arbitrarily chosen equal to each of the spheroidal 
mode resonances. One finds a maximum change of 10% in the SNR for a resonator mistun- 
ing of about 10%. It is possible to optimize the tuning by shifting the resonator resonance 
frequency as much as indicated by the maximum of the SNR in figure [IBa. This proce- 
dure can be repeated several times. The result is shown in figure [03b, where the SNR 
is maximum for each resonator around the new set of the resonators natural frequencies 
uj\ = 27r[1047, 971, 980, 970, 960, 1027]. Such an optimal frequency set derives from a combi- 
nation of multiple coupling, modes splitting and transducers position. Once the bare sphere 
spheroidal modes and the resonator masses and positions are known, one can always find an 
optimal natural resonance for each resonator. 

We remark here that, due to the presence of so many modes and the related multiple 
splitting, it could be difficult in practice to determine the resonators frequency with an 
accuracy better than 10% and one has probably to accept a not optimized detector. The 
loss in sensitivity is in any case less than 10% for a 100Hz mistuning of the mechanical 
modes of a detector resonating at 1kHz and with a tuned matching network. 

So far we consider the transducers operating on the sphere all identical except for their 
main resonance frequency. Here we evaluate the effect on the detector anisotropy in the 
sensitivity and bandwidth, which derives from using six not identical transducers to read- 
out the five quadrupolar mode of the sphere. We arbitrary modified the parameters of each 
transducer such as mass, mechanical and electrical quality factor, transducer capacitance, 
electrical coupling factor and SQUID noise of a maximum of 15% from their optimal value. 
As shown in figure [TOjb, the detector SNR, for a circularly polarized GW burst with Tqb = 
3fiK, is reduced about the 30%. This is mainly due to a decrease in signal bandwidth as 
one can see from figure [TTlb. 

In the following we study how the detector sensitivity decreases when the transducer 
electrical resonator is not perfectly matched to the mechanical resonator. Such a situation 
could arise in practice when the electrical mode cannot be arbitrarily adjusted to the opti- 
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FIG. 16: Resonators mistuning. The graph shows the detector SNR for a Tqb = S^iK one-cycle cir- 
cularly polarized GW burst as a function of the mistuning parameter uj /ojq for each resonator in the 
case of a) a starting arbitrary set of natural resonances ojq = 27r[987, 1001, 1001, 1008, 1012, 1017] 
and b) an optimized set uj\ = 27r[1047, 971, 980, 970, 960, 1027] obtained after few tuning iterations. 

mum value define by Eqs. (1381) as a consequence, for example, of a voltage leakage in the 
bias lines. One finds that, when the electrical mode of only one transducer is mistuned, 
even up to about 30% of the optimal frequency, very little effect is observed in the SNR and 
bandwidth. This is evident from the full circle data of figure [T71 where the detector SNR 
and bandwidth for a Tqb = 3/j.K circularly polarized GW burst is given as a function of the 
mistuning parameter u/u e i, being u e i equal to the optimized mechanical resonator frequency 
set derived above. When all the electrical modes are decoupled from the mechanical modes, 
the SNR and bandwidth decrease of a factor proportional to the mistuning factor as shown 
by the open squared data of figure [T71 

In the numerical analysis of resonant detectors one always considers the mechanical res- 
onators as point-like masses. A real transducer is sampling the sphere surface radial dis- 
placement in many points belonging to the contact surface between the resonator spring 
and the sphere. The actual displacement can be seen as the results of an average of such a 
sampling. We calculated here the detector sensitivity when a real resonator is considered. 
Three kinds of resonators are generally used in resonant detectors: mushroom, membrane 
and "rosette" resonators. All of them have in common the fact that the spring is attached 
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FIG. 17: Mistuning of the electrical mode. The graph shows the detector SNR (a) and bandwidth 
(b) for a Tqb = SfiK one-cycle circularly polarized GW burst as a function of the mistuning 
parameter to/ujei where ui e i is equal to the optimized mechanical resonator frequency set. The data 
show the effect when the electrical mode of only one transducer is detuned (full circles) and when 
the electrical modes of all the six transducers are equally mistuned (open squared). 

to a support, which is rigidly connected to the antenna's surface. The variable read-out 



considered in our calculation, in analogy with 34|], can be written as 



x(t) = J f(rt,r T )y(Tt,r T ,t)dr t dr r (39) 

where y(r,t) is the radial displacement of the intersection points between the support and 
the springs of the mechanical resonator used to amplify the sphere displacement, /(r) is a 
form factor and the integral is calculated along the resonator support-spring intersection. 
We consider transducers with cylindrical symmetry where r t is the radius of the support- 
spring intersection and r r gives the radial position of the intersection points referred to the 
sphere centre. For "mushroom", "rosette" and "membrane" resonators those intersection 
are, respectively, the "mushroom" leg section, the cylindrical section where the membrane is 
attached to the resonator ring support, and the "rosette" spring sections at the attachment 
point with the ring support. 

The form of the radial displacement depends strongly on the resonator geometry and 
springs topology, making it difficult to find an analytical expression like the one for the 



gaussian laser beam considered in 



341 ]. and 35|. In the case when only the first spheroidal 
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modes are considered, the form factor for a "membrane" and "rosette" transducer becomes 

fmemb (r t , T r ) = — , 40 

l-nrMr A 12 {r s ) 

where 8r r is the spring thickness and A±2(r) is the spheroidal quadrupole radial amplitude 



function described in 



111 ]. For a "membrane" transducer s — 1 , while for a "rosette" 



transducer s = S ros /S mem b < 1 is the ratio between the intersection surfaces defined by the 
"rosette" springs and a membrane of the same thickness. 
In the case of a "mushroom" transducer, we have 

_ 1 A 12 {r r ) 

J mush oa I \ ' V / 

In proximity of the sphere surface the radial amplitude is slowly changing and the ratio 
Ai2{r r ) I 'Ai 2 {r s ) can generally be approximated to 1. 

In figure [T8l the SNR and bandwidth for a circular polarized Tqb — 3{iK GW burst of a 
nearly quantum limited sphere with six resonators is shown as a function of the resonator 
radius r t . The improvement in sensitivity, obtainable as described above by using massive 
resonator, with \i = 0.01, is slightly reduced by the resonators large radius r t ~ 0.2m. The 
SNR decreases of about 10% in this case. A "mushroom" resonator is preferable, as shown in 
figure dH but high Q massive resonators are difficult to achieve for such a geometry. In the 
sky maps shown in figure [TjH obtained from 1000 randomly distributed events of circularly 
polarized GW bursts, the SNR and bandwidth anisotropy for a transducer with r t = 0.2 m 
is shown. One finds up to 15% of asymmetry in the SNR and up to 40% in the detector 
bandwidth. 



V. CALIBRATION OF A SPHERICAL GRAVITATIONAL WAVE DETECTOR 

To calibrate a gravitational wave detector one needs to postulate a model of the complete 
detector whose parameters are experimentally determined. The calibration is performed in 
two steps. First the effective temperature of the modes is estimated. This is important in 
order to understand the detector dynamics and to estimate the noise contributions. Sec- 
ondly, one has to measure the detector response to an applied force, which is equivalent to 
evaluating the transfer function G s i gj i in Eqs. ( T261) . We stress here that a full antenna char- 
acterization implies the estimation of the terms of the matrix G s i gi i, which relates the input 
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FIG. 18: Detector SNR and signal bandwidth as a function of the resonator radius rt for a 
"membrane" resonator with mass ratio [x = 0.005, open squares, and [i = 0.01, full circles, and a 
"mushroom" resonator with fi = 0.01, open triangles. 




FIG. 19: Detector SNR and signal bandwidth anisotropy for a circular polarized GW burst with 
Tqb = 3{iK when a "rosette" resonator with rt = 0.2m is used to read out the sphere quadrupolar 
modes. 

signal and noise to each transducer SQUID input current. It can be achieved by injecting 
a known signal at each step of the transducer chain, i.e: radial forces to the five spheroidal 
modes of the sphere, radial forces to each resonator, voltage signal in the superconducting 
transformer and voltage signal at the SQUID input. In order to evaluate the sensitivity of 
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the antenna to GW signal, however, one needs especially to measure the transfer function 
which converts the forces applied to the spheroidal modes to the SQUID input current of 
each transducer. By measuring the noise of each transducer during normal operation, one 
can finally estimate the strain sensitivity as in Eqs. (|26p . 

We describe below two methods that can be used to perform a complete calibration of a 
spherical detector in analogy with the method used for the bar detector AURIGA l|, 

A. Modes equivalent temperature 

According to the fluctuation dissipation theorem, the voltage and current power spectra 
observed in an electrical circuit, when only thermal noise sources are presents, are given by 

S v = AkTRe{Z(u)), S T = AkTRe (^r-A , (42) 

\Z(lo)J 

where Z[iS) is the impedance seen at the input of the amplifier used for the read-out and T is 
the thermodynamic temperature. By measuring the transducer output spectrum in normal 
operation and the input impedance of the circuit one can evaluate the equivalent temperature 
of the transducer chain. If only thermal noise is present, the equivalent temperature of the 
chain should be equal to the thermodynamic temperature of the experiment. Because the 
SQUID is not an ideal amplifier, its current and voltage noise give a contribution to the total 
transducer output noise. This contribution is in general not negligible. The total current 
noise measured at the input of each transducer line SQUID can be written as l|, 

Sj = AkT^LRe (-L-) + + Sm(T), (43) 

Qei \Z(u)J \Z(u)\ 2 

where the second term on the right hand side of the equation is the back action noise 
contribution from the SQUID amplifier with S vv given by Eqs. ( 1361) . Q a is the apparent 
quality factor produced by the damping and Q e \ is the intrinsic quality factor of the electrical 
matching network. The factor §r m the thermal noise appears when the cold-damping 

Wei 

network is active in the read-out circuit lj. It comes of the fact that the damping is only a 
lossless electronic feedback effect and there is no dissipation associated to it. From Eqs. (143]) 
we see that when the SQUID noise parameters S VV (T), and Sa(T) are known, and the 
impedance Z(cu) is measured, from the fit of the SQUID output current one can estimate as 
a fitting parameter the equivalent temperature of each transduction chain. One can measure 
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the input impedance Z(u) seen by the SQUID by injecting a sinewave signal with defined 
level through a calibration coil weakly coupled to the SQUID input circuit. 

We can simulate numerically the calibration procedure. Denoting by M ca i the mutual 
inductance between the calibration coil and the SQUID input circuit and by R ca i the resis- 
tance of the calibration line, to estimate the impedance Z{iS) one has to solve the system 
of equations in Eqs. (fTUl) . All the terms on the right hand side are zero but the SQUID 
input voltage, which is given by V n = ^^ Ical V ca u with V ca i the voltage of the injected cali- 
bration signal. By measuring the output response of the SQUID amplifier we obtain a direct 
estimation of each transducer admittance from the following weighted average: 

1 Re{I n V:} + jIm{I n V*} 



M IK 



GT, (44) 



where I n is the current at the input of each SQUID and G™ was defined in Eqs. (12 lj) . The 
weighted average gives a more precise result than simply measuring I n /V n . We remark that 
both amplitude and phase of the input and output signal have to be measured to estimate 
Z(u). 

The admittance 1/Z(u) can be approximated by a complex polynomial expansion with 
N p — 5 + 2 x N poles and ^ = 5 + 2x^-1 zeros as follows 

1 = /^^-ft)^-^ ^^) (45) 

If Uk and Q a ,k are the resonant frequency and the apparent Q-factor of each measured 
resonance oft he system, we have pt = juj k — u>k/2Q a> k-The zeros qu can be written in the 
same form and have frequency and Q which depends on the modes coupling. Once the 
admittance is measured, a polynomial fit can be performed in order to find its zero and 
poles. 

We simulated numerically the mode temperature calibration procedure. The current 
noise at the input of the SQUID amplifier was estimated from Eqs. (1421) by calculating the 
impedance seen by the SQUID in the form of the admittance 1/Z(u). For simplicity we 
only consider two resonators located in position 1 and position 2 of the TI arrangement. 
The detector is at T = 100 mK with &T/Q = 2.5 x 10 -7 , a SQUID energy resolution of 
200H and an apparent quality factor Q a = 600. The electrical modes of each transducer are 
tuned to the mechanical modes. 
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FIG. 20: Mode temperature calibration. The current noise density at the input of transducer 1 
SQUID is estimated by calculation the admittance matrix l/Z(to). The contribution from thermal 
noise b), back-action noise c) and total noise including the additive SQUID current noise a) are 
shown. The dashed curves are obtained from the pole-zero polynomial expansion. 

Figure [201 shows the SQUID input current noise contribution from the thermal and back- 
action noise and SQUID additive current noise of transducer line 1. The zero-pole approx- 
imation of Eqs. (145|) is shown with dashed curves on top of the simulated current noise 
contributions. 

B. Force calibration and transfer functions measurements 

To calibrate a spherical detector one should be able to experimentally evaluate the transfer 
functions which relate the output current with the GW force acting on the spheroidal modes. 
When a force is applied at a point (8 C , <j) c ) on the surface of the sphere, all the 5 + 2N modes 
of the sphere and transducers are excited at an amplitude which depends on the calibrator 
position, the resonators positions and the coupling between resonators and spheroidal modes. 
A calibrator, consisting of a capacitive transducer with the main resonance frequency u> ca i 
not tuned with the fundamental frequencies u of the detector, can be used to convert an 
electrical signal into a constant force acting on the sphere modes. For u ca i » uo, such a 
force is given by 




ca I 



(46) 
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where C ca i, E ca i and m ca \ are respectively the capacitance, the dc bias electric field and 
the mass of the calibrator. V ca i(u>) is the excitation voltage given to the calibrator at the 
resonance frequencies of the detector. 

A force acting on the sphere surface always excites a combination of the 5 spheroidal 
modes. If denote by F c = (F 1 ...Fj^ c ) a vector of N c radial forces applied to the N c points 
P c = (&c, 0c), with c = l...N c , the corresponding forces F m acting on the 5 spheroidal modes 
can be described by the calibrators pattern matrix Be: 

F m = aB c F c . (47) 

From Eqs. we find Bc m , c — '^m,c(^c,4 ) c)- The pattern matrix B c has the same physical 
meaning as the matrix B, but is referred to the calibrator positions. From Eqs. ( I4T1) it is 
clear that the gravitational forces acting on each of the five quadrupole modes, distributed 
over the entire volume, can be simulated by a linear combination of radial forces acting on 
the N c ^ 5 calibrators. This conclusion is based on the assumption that the real spherical 
antenna dynamics is the same as that of an ideal sphere that can fully be modelled using 
the spherical harmonics approach described above. The validity of this assumption could 
be evaluated through a finite element analysis (FEA) of the detector structural model or 
by experimentally measuring the six transducers response to several linear combinations of 
excitation given through the calibrators. 

To operate the sphere as a GW detector, only one force calibrator, mounted at an arbi- 
trary position on the sphere surface, is necessary when N = 6 transducers are used in the 
TI configuration. This is due to the existing one-to-one relation between the mode channels 
and the forces acting on the spheroidal modes. It is convenient to proceed by transforming 
the N measured transducers outputs into the 5 mode channels, which directly describe the 
GW amplitude acting on the spheroidal modes In such a framework we derive the 

optimal filter for each mode channel to define an operative procedure to signal detection 
with a GW spherical resonant antenna. In order to do this we define the 5x5 mode channel 
noise spectral density matrix as 

Si, B (w) = BS I (o;)Bt, (48) 

and, using the admittance matrix G s i gi i introduced in Eqs. (1261) . the 5x1 mode channel 
signal vector 

I g)Sig = Bli = aBG sig iB c F c . (49) 
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By applying a known constant force to the calibrator F c , using Eqs. (147|) we get the linear 
combination of spheroidal modes forces generated by the calibration force. By measuring 
simultaneously amplitude and phase of the six transducer outputs we obtain the mode chan- 
nel response from Eqs. fT4^]) . and, using Eqs. (1261) . estimate the detector strain sensitivity. 
The matrix Si^u;) in Eqs. (I4"gj) is diagonal because the mode channels are statistically 
independent. The only assumption made here is that the the pattern matrix B derived 
above for an ideal sphere can be used to describe the dynamics of the real detector with 
six transducer in the TI arrangement. This can be experimentally verified in a separate 
cryogenic experiment, using for example a set of six calibrators in the TI configuration with 
a 60° de-phase of the azimuthal angle with respect to the transducer arrangement. Each 
spheroidal mode force can be reproduced by simultaneously applying a linear combination 
of constant forces on the six calibrators, as described by Eqs. (T4"Tj) . When N < 5 transducers 
are used we cannot transform the transducer outputs in mode channels. A set of at least 
5 calibrators is then necessary to experimentally measure the transducer response to each 
mode channel. If only one calibrator is used, the detector output can only be calibrated for 
a particular combination of the five quadrupolar forces F m . 

Each element on the diagonal of Sj™^ = I m (u}I^(oS) can be written in terms of polynomial 
ratio where the poles are the same as the ones derived in the impedance measurement of 
each transducers, see Eqs. ( |4"5]1 . and the zeros depend on the current function I m (u). From 
the factorization of I m (u) one finds 

U«) = ^n M {ju ,_ pk){ju ,_ pl) . (50) 

where So, m is the wideband noise of the m th mode channel, and equals the amplifier additive 
white noise if all the transducers SQUIDs are identical. From the definition in Eqs. (1501) . 
S}^ is real for real oj and the number of zeros and poles is the same as a consequence of 
assuming S^to be purely white. 

The transfer functions for a GW signal, which convert the quadrupolar modes forces 
into mode channels currents according to Eqs. (j4"9l) . contain the same poles {pk} and their 
factorization becomes 

H m {uj) = H nhcal {u) — — . (51) 

In the equation above N p > N r and H mjCa i(uj) is a force calibration constant which has to 

42 



be experimentally determined at each cool down. 

From now on we apply to each mode channel the standard Wiener-Kolmogorov (WK) 
filtering operation, developed so far for bar detectors [ijj]. We follow here the approach de- 



scribed in 23j. The best linear estimate of the amplitude ho of a given signal h(t), with max 
h(t) = ho at the arrival time t = 0, buried into an additive, zero mean, stationary gaussian 
noise 77 can be obtained by correlating the mode channel output Y m (u>) = h(u)H m (u) +i]{ui) 



to the matched WK filter 



39|, 



W m (u) = a A , (52) 

where a\ = J duj\H m (uj)h(uj)\ 2 / Sj 17 ^ (00) is the variance of the noise after the filtering. The 
WK filter splits as the product L m {uj)M m (uj)h*{uj), where L m (cu) is a whitening filter for 
the noise S/"g(u;) given by 



10J 



T-l (l ,\ _ 0-l/2jrN P (jv-Pk,m)(3U-P*k,m) 

(juj - qk, m ){ju - qu m ) 



M m is a bandpass filter around the frequencies = |-^w(?fc,m)| and bandwidths Auk = 
2Re(qk, m )- Such a bandwidth is generally much larger than the intrinsic bandwidth of each 
resonance and, in the case of a transducer with coupled electrical modes, it can reach values 
as large as ~ 200Hz. One finds 

M m {u) = a A ■ S Q>r i — (54) 

Rk^iiju + qk,m){jv + qlJ 

The product L m (uj)M m (u!) is the WK filter for a delta-like GW pulse. The extra term h(u) 
should be added when a general GW signal h(t) is considered. It can be shown that for a 
resonant bar detector, N r = 1 and n = 0. In this case the WK filtering procedure is then 
fully defined by the zeros and poles p k of the noise power spectrum Sj(uj), the additive 
amplifier white noise So and the calibration constant H (uj). For a spherical detector this is 
true only if the quadrupolar modes degenerate into a single resonant frequency and all the 
transducers have the same resonance. When multiple resonances are present, as is the case 
for real spherical detectors, the zeros in Eqs. ( 15TT) do not cancel i.e., N r > 1 and rk, m 7^ 0. 
This is due to the fact that when the quadrupolar modes are non-degenerate, a mixing 
occurs between mode channels. A fraction of the signal which should only go to one mode 
channel leaks into the others. The WK filtering procedure should include the extra set of 
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FIG. 21: Modes channels response from a simulated linearly polarized burst coming along the z 
direction in the detector frame for a detector with degenerate, a), and non-degenerate, b), modes. 
For a detector with degenerate modes, a), only the second mode gives a larger signal than the 
noise. For a real detector with non-degenerate modesb, a mixing between the modes is present 
around the resonances making the analysis more complex 

parameters rf. jm , whose total number depends on the mode and transducer considered and 

must be experimentally determined. 

In figure [21] the mode channels response is plotted for a linearly polarized burst coming 

along the detector z direction. For a detector with degenerate modes, figure [2TJa, only the 

second mode gives a larger signal than the noise, as can be expected for a signal from that 

particular direction and polarization. For a realistic detector with non-degenerate spheroidal 

modes a mixing between the modes is present around the resonances making the analysis 

more complex, see figure [2Tlb). However, the energy stored in the second mode, which can 

be derived from the integral of the SNR density, is larger than in the others and the incoming 

wave direction can still be reconstructed without any significant accuracy loss with respect 

to the degenerate case. To estimate the incoming wave direction we can use the approach 
□ 

derived in [12J, using standard theory of signal detection. After being optimally filtered, the 
5 mode channels generate a set of 5 amplitudes, g m , for an input GW burst. The likelihood 
function for a detector with stationary and Gaussian noise is given by 

where h m is the expected gravitational wave signal amplitude, g m is the mode channels 
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FIG. 22: Signal direction reconstruction. An overall sky search is performed calculating the like- 
lihood function A. We applied simulated burst signals with SNR ~ 30, 300, 1000, amplitude 
/io = 0.34 • 10 -20 , 1.1 • 10~ 20 , 2 • 10 -20 and tqb ~ lms, linearly polarized and coming from a source 
at declination 5 = 20° and ascension a = 70°. A single spherical detector cannot distinguish 
between sources laying in two opposite hemispheres. 

amplitude obtained after the WK filtering procedure described above and a m is the variance 
of g m . A likelihood map is generated by plotting this function in the decimation- ascension 
plane, (5, a). The maximum value of A gives to the estimated wave direction. In the example 
below, a realistic detector at the quantum limit, with optimized parameter is considered. 
We applied simulated burst signals, 1ms long and with amplitude ho = 0.34 • 10 -20 , 1.1 ■ 
10~ 20 , 2 • 10~ 20 and SNR ~ 30, 300, 1000, with linear polarization h+ = h and h x = 0, 
coming from a source at declination 5 = 20° and ascension a = 70°. Figure [22] plots the 
resulting likelihood functions. 

The accuracy in the direction and polarization estimates depends on the SNR and is 
equal to Af2 = 2-k/SNR [3], [l3]; it is direction independent when a sphere with more 
than 5 identical, point-like transducers is considered 44j. As it comes clear from the plot, 
a single spherical detector cannot distinguish between sources laying in the two opposite 
hemispheres. In the example described above we assumed for simplicity that the signal 
polarization were known. When this is not the case, one should include also the first Euler 
angle ip as a variable to estimate the likelihood ratio. 

As discussed in [13] , for a known signal, the optimal detection strategy for a vector output 
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detector like a sphere with N transducers is to compute the value of the optimal linear filter 
with scalar output v which maximizes the SNR. This procedure is useful if one wants to 
monitor the total energy in the sphere and claim a detection only on the base of an excess 
of absorbed energy. No direction information is possible in this way. 

If we let I = I\...In be the transducers output stream, a linear filtering operation is 
performed which, in frequency domain, can be described by 

v (u) = W + (w)I(u;), (56) 



where W is a vector transfer function which maximizes the SNR and was found 13J, to be 



WH=4HS, 1 H. (57) 



In the above, I s i g = G S i g (o;)F^ 1 is the current generated in each transducer by a force 
on the sphere generated by a GW signal. As shown in 12|, the maximum SNR for a 
multichannel detector is the sum of the maximum SNR of each individual channel. The 
optimal linear filter introduced above for each mode channel is then useful in order to 
compute the scalar output v which maximizes the SNR. 



VI. CONCLUSIONS 



We derived a complete and detailed electro-mechanical model for a spherical gravitational 
wave detector operating with multiple, two-mode capacitive transducers where the electrical 
resonant mode of a superconducting matching LC resonator can be tuned to the resonant 
modes. The signal current form the matching network is read out by SQUID amplifiers. 
The model allows to numerically calculate the sensitivity of a realistic detector and to study 
in detail the effects of the main mechanical and electrical parameters of the displacement 
read-out system on the strain sensitivity and bandwidth. All the known noise sources are 
discussed and considered in the model. 

A complete numerical analysis has been performed for a 2 meter in diameter, 30 ton 
in mass, CuAl spherical detector operating at ultracryogenic temperatures. We derived 
the sensitivity for the spherical detector in its initial phase of development, when a single 
transducer is used, and when the detector operates with six transducers and becomes fully 
omnidirectional. The sensitivity is evaluated when the detector operates by making use of 
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available technology and when it works at the quantum limit. We have shown that, in order 
to improve the strain sensitivity towards the quantum limit one should operate the detector 
at temperature of about T = 20mK with electrical and mechanical quality factor as high as 
5 ■ 10 7 and massive mechanical resonators. 

Direction anisotropics in the detector sensitivity and signal bandwidth are studied for a 
not-ideal detector operating with not identical, partially tuned and real-size resonators. The 
models made so far always consider rather generic, an point-like transducers, neglecting the 
fact that those sensors are in practice rather large and spatially distributed on a significant 
fraction of the sphere surface. We investigated here the validity of such an assumption. 

Finally we described and numerically verify a complete calibration procedure, which 
makes use of techniques available for bar detectors. Similar algorithms can be used on 
spherical detector for diagnostic purpose and to derive the direction and polarization infor- 
mation from the detected signal. 
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